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Abstract  Matched  field  processing  (MFP)  has  been  extensively  explored  for  use 
in  detecting  and  localizing  underwater  sources.  By  matching  the  received  data 
to  replica  vectors,  which  account  for  multipath  effects,  three-dimensional  source 
localization  is  possible.  A  significant  problem  for  MFP  is  that  the  beamformer 
output  often  displays  high  ambiguities  in  range  and  depth,  with  a  single  source 
leading  to  multiple  spatial  peaks.  This  is  particularly  true  for  arrays  with  limited 
vertical  aperture  because  of  their  limited  ability  to  resolve  multipath.  Many  of  the 
ambiguous  regions  may  be  essentially  in  the  main  lobe  of  the  cell  of  interest, 
making  adaptive  suppression  of  them  difficult. 

In  this  paper  we  explore  methods  for  analyzing  and  exploiting  the  structure  of 
the  ambiguity  surfaces.  Statistical  clustering  methods  allow  regions  of  similar 
response  in  the  ambiguity  surface  to  be  identified,  using  the  generalized  cosine 
between  replicas  as  a  measure  of  distance.  The  cluster  analysis  gives  a  map 
of  ambiguous  regions  in  the  MFP  search  space.  We  present  several  ways  of 
exploiting  this  map.  First,  we  demonstrate  that  computational  gains  are  possible 
by  treating  each  cluster  as  a  separate  "beam"  that  is  represented  by  a  single 
replica  vector.  Because  the  ability  to  resolve  sources  is  a  function  of  SNR,  the 
appropriate  cluster  size  will  depend  on  SNR  as  well.  This  can  be  built  into  the 
adaptive  beamforming  by  allowing  beamforming  at  multiple  resolutions.  Second, 
we  show  how  the  cluster  map  can  be  used  in  post-processing.  By  "peak-picking" 
in  cluster  space,  rather  than  in  spatial  dimensions,  we  are  able  to  identify  and 
discard  spatial  ambiguities  that  result  from  the  presence  of  a  strong  source. 
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ABSTRACT 

Cluster  analysis  provides  a  tool  for  mapping  out  regions 
of  ambiguous  response  in  sparse  array  beamforming  prob¬ 
lems.  This  paper  discusses  clustering  and  its  application  to 
matched  field  processing  (MFP)  problems  in  ocean  acous¬ 
tics.  The  map  of  the  ambiguity  volume  provided  by  clus¬ 
tering  can  be  used  for  improved  interpretation  and  post¬ 
processing.  By  peak-picking  in  cluster  space,  rather  than 
in  spatial  dimensions,  we  are  able  to  identify  and  discard 
ambiguous  peaks  that  result  from  the  presence  of  a  strong 
source.  As  shown  in  a  data  example,  clustering  can  help 
collapse  3-D  MFP  output  to  bearings-only  while  preserving 
signal  gains  obtained  by  accounting  for  multipath.  Cluster¬ 
ing  can  also  provide  computational  gains  by  allowing  elim¬ 
ination  of  highly  redundant  replicas. 

1.  OVERVIEW  OF  MATCHED  FIELD 

Matched  field  processing  has  been  extensively  explored  for 
use  in  detecting  and  localizing  underwater  sources.  MFP 
uses  ocean  propagation  models  to  account  for  multipath  when 
generating  replica  vectors.  If  accurate  environmental  in¬ 
formation  is  available,  sources  can  be  localized  in  range, 
depth,  and  bearing.  Accurately  modeling  multipath  can  also 
yield  mismatch  reduction  and  detection  gains  as  compared 
to  direct-path  beamformers.  A  significant  problem  for  MFP 
is  that  the  beamformer  output  often  displays  high  ambigu¬ 
ities  in  range  and  depth,  with  a  single  source  leading  to 
multiple  spatial  peaks.  This  is  particularly  true  for  arrays 
without  significant  vertical  aperture  because  of  their  limited 
ability  to  resolve  multipath.  MFP  for  these  arrays  can  be 
thought  of  as  a  form  of  sparse  array  processing  with  resul¬ 
tant  high  ambiguities.  Many  of  the  ambiguous  peaks  may 
be  within  a  few  dB  of  the  overall  peak,  so  they  are  not  easily 
suppressed  with  adaptive  beamforming. 

This  work  was  sponsored  by  DARPA-ATO  under  Air  Force  Contract 
F19628-00-C-0002.  Opinions,  interpretations,  conclusions,  and  recom¬ 
mendations  are  those  of  the  authors  and  are  not  necessarily  endorsed  by 
the  United  States  Government. 


MFP  replicas  are  generated  for  a  set  of  look  directions 
in  range,  depth,  and  azimuth.  If  we  define  0  =  (r,  0,  z)  as 
the  look  direction,  the  MFP  output  for  an  N  element  array 
at  frequency  /o,  time  to,  and  direction  0  is  given  by 

P(fo,to,@)=  \wH(f0,t0,<d)x(fo,to)\2  (1) 

where  w(/o,  to,  0)  is  the  N  x  1  adaptive  weight  vector  and 
x^(/o,  to)  is  the  current  data  snapshot.  Applying  the  adap¬ 
tive  weights  to  individual  snapshots,  rather  than  to  a  covari¬ 
ance  matrix,  helps  reduce  smearing  due  to  target  motion. 
The  adaptive  weight  is  based  on  the  corresponding  replica 
vector  v(/o,  to,  0)  and  the  N  x  N  sample  covariance  ma¬ 
trix.  The  replica  vectors  are  calculated  using  an  appropriate 
propagation  model  [4] ;  an  adiabatic  normal  modes  approach 
is  used  here.  The  sample  covariance  matrix  is  computed  us¬ 
ing  a  time  average  of  FFT  snapshots.  In  passive  sonar  the 
covariance  matrix  contains  data  from  the  target  as  well  as 
interferers,  so  signal  mismatch  and  self-nulling  can  be  sig¬ 
nificant  problems. 

In  this  paper  we  will  apply  the  dominant  mode  rejection 
algorithm  with  white  noise  gain  control  (DMR-WNGC)  for 
adaptive  processing  [2].  This  algorithm  approximates  the 
sample  covariance  matrix  as  the  sum  of  a  dominant  eigen¬ 
vector  subspace  and  a  noise  subspace.  Protection  against 
mismatch  is  provided  by  adding  a  diagonal  load  that  varies 
with  look  direction  to  satisfy  a  white  noise  gain  constraint. 


2.  CLUSTERING  CONCEPTS  AND  APPROACHES 


Statistical  clustering  methods  [3]  enable  identification  of  re¬ 
gions  of  similar  response  on  the  ambiguity  surface.  Cluster¬ 
ing  algorithms  require  a  distance  measure  to  estimate  the 
similarity  between  two  given  vectors.  For  our  problem,  a 
natural  distance  measure  is  given  by  the  cosine  between  two 
replicas  [1].  The  distance  between  two  frequency  domain 
replica  vectors  and  vy  is  given  by: 


dist(vj,  Vj)  =  1  - 


2  ' 


(2) 


This  distance  metric  is  then  used  to  determine  which  vectors 
should  be  clustered  together.  Clustering  has  been  previously 
applied  to  sparse  array  direction-finding  problems  using  the 
same  distance  metric. 

Algorithms  used  to  perform  clustering  have  two  main 
outputs.  The  first  is  a  set  of  cluster  centers ,  each  of  which 
is  a  single  replica  vector  that  is  chosen  to  be  representa¬ 
tive  of  all  members  of  its  cluster.  As  shown  in  Figure  1, 
the  cluster  center  can  be  considered  to  be  equivalent  to  the 
axis  of  main  response  (MRA)  for  a  standard  beampattern. 
The  second  output  is  a  cluster  map  listing  all  of  the  repli¬ 
cas  which  belong  to  the  cluster  represented  by  each  center. 
These  points  correspond  to  main  lobe  look  directions.  For 
the  methods  used  here,  the  user  must  specify  a  cluster  ra¬ 
dius,  or  maximum  allowable  mismatch  between  each  clus¬ 
ter  center  and  all  cluster  members.  This  radius  is  equivalent 
to  the  “dB  down”  number  typically  used  to  define  main  lobe 
width.  The  utility  of  clustering  is  indicated  by  the  sparse  ar¬ 
ray  beampattern  shown  in  Figure  1 .  When  grating  lobes  are 
present,  look  directions  that  are  within  the  specified  clus¬ 
ter  radius  may  be  spatially  distributed  over  a  wide  region. 
Clustering  provides  an  automated  way  to  identify  these  dis¬ 
tributed  regions  of  similar  response. 


Figure  1:  Beampatterns  for  nominal  well-sampled  (solid 
line)  vs.  sparse  arrays  (dashed  line). 

Commonly  used  clustering  algorithms  such  as  k-means 
or  heirarchical  clustering  are  extremely  slow  when  applied 
to  large-dimensional  problems  such  as  MFP.  We  have  found 
that  the  following  ideas  help  significantly  for  the  MFP  clus¬ 
tering  problem: 

•  the  clustering  method  should  be  “greedy,”  meaning 
that  it  does  not  seek  the  optimal  clustering  solution; 

•  a  divide-and-conquer  (denoted  DnC)  approach  should 
be  used,  in  which  the  method  operates  on  sub-sets  of 
replica  vectors  to  build  up  the  solution; 


•  if  the  array  geometry  and  environment  vary  smoothly 
with  time,  techniques  for  quickly  updating  previous 
cluster  results  can  be  employed. 

In  test  cases,  using  a  greedy  algorithm  gave  a  speedup  of 
roughly  lOOOx .  Combining  the  DnC  idea  with  a  greedy  al¬ 
gorithm  gave  another  10-20  x  speedup,  while  fast  updating 
techniques  gave  an  additional  speedup  of  3-5  x . 

3.  CLUSTER  ANALYSIS  OF  MFP  REPLICAS 

In  this  section  we  will  use  cluster  analysis  to  gain  insight 
into  the  3-D  structure  of  the  MFP  array  manifold.  MFP  sim¬ 
ulations  were  carried  out  using  environmental  parameters 
from  the  Santa  Barbara  Channel  (SBCX)  experiment  [6]. 
Simulations  were  done  comparing  a  notional  single-line  ar¬ 
ray  (101  phones  at  2  m  spacing)  to  a  volumetric  array  con¬ 
sisting  of  three  vertically  stacked  arrays  with  30  m  total  ver¬ 
tical  aperture.  Simulations  were  done  for  frequencies  and 
source  levels  used  in  the  actual  SBCX  test  (12  tones  from 
64-338  Hz  at  159  dB  re  1  /iPa  )  with  80  dB  white  noise 
added. 

Figure  2  shows  how  the  members  of  an  example  clus¬ 
ter  at  94  Hz  are  distributed  in  three  dimensions.  Clusters 
were  generated  with  1.5  dB  cluster  radii  for  the  single-line 
array  described  above.  Cluster  members  are  well  grouped  in 
bearing  but  are  widely  distributed  in  range  and  depth.  While 
not  shown  here,  similar  patterns  are  seen  when  plane- wave 
replicas  are  used  as  cluster  centers  for  clustering  the  MFP 
replicas.  When  vertical  aperture  is  added  the  cluster  mem¬ 
bers  become  more  clumped  in  range  and  depth.  This  plot 
demonstrates  several  important  features: 

•  When  vertical  aperture  is  limited,  source  azimuth  is 
much  better  resolved  than  range  or  depth.  Thus  the 
SNR  and  environmental  calibration  requirements  for 
correctly  localizing  source  bearing  will  be  much  lower 
than  for  range  /  depth  localization. 

•  While  cluster  members  for  the  single-line  array  are 
widely  scattered  across  range  and  depth,  several  clus¬ 
ters  will  be  required  to  fully  sample  the  search  space 
along  a  particular  bearing. 

To  further  explore  the  properties  of  MFP  clusters  we 
consider  an  endfire  search  problem  for  the  single-line  array. 
The  array’s  sensitivity  to  multipath  is  highest  at  endfire,  so 
the  potential  benefit  of  accounting  for  multipath  propaga¬ 
tion  is  greatest  in  this  region.  As  a  proposed  search,  we 
consider  a  volume  extending  from  -30  to  +30°  around  for¬ 
ward  endfire,  over  the  full  water  depth  of  0-210  m,  and  from 
1-10  km  in  range.  The  space  is  gridded  with  15  azimuthal 
beams  and  is  sampled  every  10  m  in  depth  and  every  50  m 
in  range. 


Figure  2:  Three-dimensional  view  of  example  cluster  for 
single-line  array,  94  Hz  simulation. 

We  define  a  coverage  metric  for  comparing  the  use  of 
direct-path  and  clustered  MFP  replicas.  Coverage  is  defined 
as  the  fraction  of  the  search  volume  vs  for  which  the  mis¬ 
match,  minimized  over  the  set  of  replica  vectors  vq  ,  is  be¬ 
low  some  threshold  7  for  allowable  mismatch: 

mindist(vs,  vg)  <  7.  (3) 

0 

High  coverage  at  a  low  7  indicates  that  aggressive  adapta¬ 
tion  is  possible  without  significant  target  self-nulling. 

To  add  realism,  environmental  mismatch  is  introduced 
into  the  MFP  replicas;  the  starting  compressional  sound- 
speed  in  the  seabed  has  been  lowered  from  1628  to  1600 
m/s,  and  sound  speed  profile  points  in  the  water  column 
have  been  perturbed  randomly  by  ±1  m/s.  This  environ¬ 
mental  mismatch  gives  roughly  1  dB  beamformer  mismatch 
across  much  of  the  search  space. 


100  150  200  250  300 

Frequency,  Hz 


Figure  3:  Percent  coverage  for  endfire  search  volume,  al¬ 
lowing  for  a  maximum  mismatch  of  7  =  1.5  dB. 

Figure  3  compares  coverage  values  calculated  for  plane- 


wave  and  clustered  MFP  replicas,  allowing  1 .5  dB  mismatch. 
The  plane-wave  replica  set  coverage  drops  from  roughly 
70%  to  20%  as  frequency  increases.  When  clustered  MFP 
replicas  are  used,  better  than  90%  coverage  is  achieved. 
Figure  4  compares  the  number  of  replicas  used  for  unclus¬ 
tered  vs.  clustered  MFP.  The  number  of  clusters  grows  with 
frequency  as  array  resolution  improves,  but  a  reduction  of 
roughly  two  orders  of  magnitude  is  seen.  Thus  the  clustered 
replica  set  achieves  increased  coverage  by  including  multi- 
path  but  requires  many  fewer  replicas  than  standard  MFP. 


Frequency,  Hz 


Figure  4:  Number  of  replicas  required  for  endfire  search. 


4.  COMPUTATIONAL  GAINS 

As  shown  above,  a  significant  reduction  in  the  number  of 
steering  vectors  can  be  achieved  by  replacing  the  MFP  repli¬ 
cas  with  the  set  of  cluster  centers.  Beamformer  output  in 
3-D  can  be  recovered  by  mapping  the  output  power  for  each 
cluster  center  to  all  cells  belonging  to  the  cluster.  This  ap¬ 
proach  achieves  computational  savings  proportional  to  the 
reduction  in  the  number  of  replicas. 

The  cost  of  this  computational  gain  is  a  potential  squint 
loss  equal  to  the  cluster  radius.  Thus  the  acceptable  squint 
loss  can  be  a  useful  criterion  for  choosing  cluster  size.  For 
the  three-line  array  simulations  described  above,  a  cluster 
radius  of  1.5  dB  did  not  introduce  significant  squint  loss  in 
reconstructed  3-D  output.  A  smaller  cluster  radius  (roughly 
0.5  dB)  was  required  to  reconstruct  single-line  results  with 
the  same  degree  of  accuracy.  Squint  loss  can  be  reduced  by 
using  ’multi-level’  clustering.  In  this  concept,  a  relatively 
coarse  clustering  can  be  used  for  an  initial  gridding  of  space. 
If  high  output  SNR  is  seen  for  a  particular  cluster  (indicating 
that  the  cluster  contains  a  source),  that  cluster  can  be  re- 
beamformed  using  either  the  actual  MFP  replicas  or  a  set  of 
smaller  clusters  generated  with  a  smaller  radius. 

The  actual  computational  gain  depends  on  the  problem 
considered.  For  the  SWellEx  data  example  shown  below 


(which  is  a  static  array  geometry)  an  overall  speedup  of 
roughly  17  x  was  observed.  When  the  array  shape  deforms 
over  time  it  is  necessary  to  re-generate  and  re-cluster  the 
replicas.  In  these  cases  clustering  can  often  reduce  the  adap¬ 
tive  weight  calculation  time  by  an  order  of  magnitude,  though 
the  overall  speedup  may  be  limited  to  3-4  x  due  to  the  need 
to  frequently  recalculate  the  replicas. 

5.  BEARINGS-ONLY  CLUSTER  LOCALIZATION 

Figure  2  above  showed  that  horizontal  arrays  typically  have 
poor  ability  to  estimate  range  and  depth  but  good  ability  to 
estimate  bearing.  The  results  also  showed  that  along  a  given 
bearing,  a  single  replica  (whether  derived  from  direct-path 
or  clustered  MFP  replicas)  did  not  provide  complete  cov¬ 
erage.  Motivated  by  these  observations,  we  will  explore 
the  use  of  clusters  for  bearings-only  localization.  The  ex¬ 
tra  information  provided  by  clustering  can  aid  in  collaps¬ 
ing  three-dimensional  search  results  to  bearings-only,  giv¬ 
ing  gains  in  track  SNR  and  deflection  ratio. 

In  general,  multiple  clusters  will  be  centered  on  each  az¬ 
imuth.  The  most  straightforward  way  to  collapse  them  is  by 
a  simple  ’max’  operation  for  all  clusters  centered  on  each 
azimuth.  Similarly,  standard  MFP  results  can  be  collapsed 
by  maximizing  along  range  and  depth.  These  maximiza¬ 
tions  will  preserve  any  signal  gain  achieved  from  accurately 
modeling  multipath.  They  will  also  increase  the  displayed 
noise  level,  an  effect  often  referred  to  as  ORing  loss. 

The  ambiguity  surface  map  provided  by  clustering  gives 
information  that  can  be  used  to  reduce  ORing  loss.  One  way 
to  separate  sources  from  noise  is  to  search  for  clusters  that 
are  peaks  in  their  local  neighborhood.  Mathematically,  we 
seek  clusters  whose  output  power  Pi  satisfies 

Pi  >  P3  Wj  e  Ni  (4) 

where  Ni  is  the  set  of  cluster  centers  j  in  the  neighborhood 
of  the  cluster  i.  While  noise  will  occasionally  give  local 
peaks  in  output  power,  the  majority  of  noise  outputs  will 
not.  Collapsing  to  bearing  by  ORing  over  the  peak  clusters 
on  each  bearing  should  therefore  give  an  improvement  in 
displayed  signal-to-noise  ratio.  The  local  neighborhood  can 
be  defined  as  the  set  of  all  clusters  whose  centers  are  mis¬ 
matched  from  the  current  cluster  center  by  less  than  some 
threshold: 

dist(v*en,v^en)  <  rj  (5) 

Setting  the  threshold  r]  to  be  twice  the  radius  of  the  individ¬ 
ual  clusters  is  a  reasonable  choice.  Several  modifications  to 
this  definition  were  found  to  be  helpful.  A  criterion  can  be 
imposed  that  the  difference  in  bearing  between  the  test  clus¬ 
ter  center  and  its  neighbors  should  not  exceed  1-2  beams. 
It  can  also  be  helpful  to  avoid  rejecting  nearby  clusters  that 
have  output  power  very  close  to  that  of  the  local  peak.  These 


modifications  help  avoid  rejection  of  a  quieter  source  near  a 
loud  source. 

This  type  of  peak-picking  operation  clearly  does  not  rep¬ 
resent  an  optimal  solution  to  the  problem  of  detecting  a  sig¬ 
nal  in  cluster  output.  The  goal  of  the  current  work  is  simply 
to  indicate  that  cluster  outputs  can  be  reduced  to  a  bearings- 
only  display  while  retaining  measurable  gains  over  a  direct- 
path  beamformer. 

Bearings-only  cluster  localization  has  been  applied  to 
data  collected  as  part  of  the  SWellEx-96  ocean  acoustics 
experiment  (www.mpl.ucsd.edu/swellex96).  For  data  anal¬ 
ysis  we  concentrate  on  event  S5,  in  which  the  towship  RJV 
Sproul  towed  a  pair  of  sources  which  each  projected  a  comb 
sequence  of  tones.  Results  are  shown  for  the  bottom-mounted 
HLA  North  array,  an  irregularly- spaced  array  of  roughly 
240  m  overall  length.  During  the  time  period  analyzed,  the 
towed  source  was  in  the  aft  endfire  (south)  direction  from 
the  array  and  was  moving  from  roughly  2  to  4.5  km  in  range. 

Figure  5  shows  a  bearing-time  record  (BTR)  for  adap¬ 
tive  processing  of  a  projected  235  Hz  tone  using  range- 
focused  replicas.  The  result  was  constructed  by  using  five 
range-focused  replica  vectors  with  focus  ranges  uniformly 
spaced  between  2-10  km.  This  represents  an  oversampling 
in  range  and  was  done  to  minimize  the  potential  mismatch 
(though  the  range-focusing  had  little  effect  near  endfire,  as 
expected).  The  DMR  beamformer  used  15  snapshots  to  es¬ 
timate  5  degrees  of  freedom,  with  a  6  dB  white  noise  gain 
constraint  applied.  While  the  high  SNR  towed  source  is 
clearly  visible,  some  time  periods  of  high  adaptive  self¬ 
nulling  are  seen  along  the  track.  At  these  time  periods  the 
signal  model  mismatch  is  high  enough  that  the  adaptive  pro¬ 
cessor  partially  cancels  the  signal  of  interest. 

Figure  6  shows  cluster  beamformer  output  where  the 
two- stage  (peak-pick,  then  OR)  method  outlined  above  has 
been  applied.  In  this  case  peak-picking  was  only  allowed 
within  ±2  azimuthal  beams  of  each  cluster  center.  The 
same  adaptive  beamforming  parameters  were  used  as  in  the 
range-focused  case.  The  cluster  result  shows  additional  power 
on  source  tracks,  due  to  mismatch  reduction,  and  a  noise 
floor  similar  to  that  seen  in  the  range-focused  result.  The 
result  is  a  display  which  has  increased  SNR  on  the  projected 
tone  and  provides  clearer  tracks  for  several  other  traces. 

The  track  detectability  can  be  measured  using  the  de¬ 
flection  ratio,  defined  as  the  difference  in  expected  values 
of  signal  and  noise,  normalized  by  the  standard  deviation  of 
the  noise: 

'/--£(Pig,}-g{Pig0>.  (6) 

a{P  |  Ho} 

This  quantity  was  calculated  along  the  track  of  the  towed 
source  for  the  cases  above.  The  expectation  was  calcu¬ 
lated  over  5  snapshots  centered  on  the  current  time  with  a 
split  window  in  bearing.  The  track  beam  and  three  guard 
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Figure  5:  SWellEx-96  bearing-time  record  (BTR),  235  Hz 
projected  tone,  for  time  period  with  R/V  Sproul  in  aft  end- 
fire.  Azimuthal  beams  are  cosine-spaced  from  —60°  to  60°. 
Results  are  shown  for  range-focused  replicas.  The  towed 
source  is  near  endfire,  in  beams  22-25. 


Figure  6:  Cluster-based  bearing-time  record  (BTR)  for  the 
event  S5  example  above.  Note  the  increased  SNR  on  the 
towed  source  as  well  as  the  improved  definition  of  the 
weaker  track  near  beam  16. 


beams  on  either  side  were  excluded  from  the  noise  calcu¬ 
lation,  with  noise  estimated  from  a  five-beam  window  on 
either  side.  To  allow  for  track  jitter,  the  deflection  ratio  was 
maximized  over  the  track  beam  and  one  beam  on  either  side. 

Figure  7  shows  the  distributions  of  along-track  deflec¬ 
tion  ratios  for  the  BTR  displays  above.  Deflection  is  also 
shown  for  a  BTR  of  MFP  results  collapsed  along  range  and 
depth  (not  shown  here).  Because  the  MFP  and  clustered 
results  avoid  significant  self-nulling,  the  probability  of  low 
deflection  is  reduced  for  these  beamformers.  The  collapsed 
MFP  result  shows  an  increase  in  overall  deflection  of  only 
1  dB  because  the  increase  in  signal  level  is  offset  by  an 
increase  in  displayed  noise  due  to  ORing  loss.  The  clus¬ 
tered  result,  which  preserves  the  signal  gains  of  MFP  but 
avoids  large  increases  in  displayed  noise,  shows  an  increase 
of  roughly  4  dB  in  median  deflection  ratio. 


Figure  7 :  Estimated  probability  distributions  for  deflection 
ratio  along  source  track,  235  Hz  tone  example  above. 


6.  USE  OF  CLUSTERS  IN  POST  PROCESSING  FOR 
3-D  LOCALIZATION 

A  final  application  of  clusters  is  in  identifying  and  rejecting 
ambiguous  peaks  in  3-D  beamformer  output.  Ideally  these 
ambiguities  would  be  suppressed  by  adaptive  processing. 
However,  the  need  to  scale  back  adaptivity  to  reduce  sen¬ 
sitivity  to  mismatch,  combined  with  the  inherent  high  side- 
lobes  associated  with  MFP,  means  that  significant  ambigui¬ 
ties  generally  remain  in  beamformer  output. 

The  cluster  map  derived  from  MFP  replica  vectors  gives 
a  way  of  representing  the  correlations  between  different  beam- 
former  steering  vectors.  By  knowing  how  the  cluster  centers 
are  related,  we  can  predict  that  a  source  in  cluster  A  will  also 
excite  nearby  clusters  B  and  C  at  some  lower  level.  If  in  the 
data  we  observe  that  cluster  A  is  more  strongly  excited  than 
B  or  C,  we  can  conclude  that  the  second  two  clusters  are 
measuring  sidelobe  energy  that  should  be  discarded.  The 
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mask  generated  by  the  peak-picking  operation  can  then  be 
used  to  zero  out  beamformer  output  for  non-peak  locations. 
This  mask  is  applied  to  the  SNR  estimate  generated  using  a 
spatial  normalization  algorithm. 

A  similar  approach  has  been  taken  in  applying  the  CLEAN 
algorithm  to  adaptive  matched  field  results  [5].  One  advan¬ 
tage  of  the  cluster  peak-picking  is  that  it  does  not  assume 
the  highest  peak  to  be  the  true  source  location,  which  may 
be  helpful  in  low  SNR  situations  when  noise  may  elevate  a 
sidelobe  to  be  the  global  maximum. 

Figure  8  shows  simulated  adaptive  MFP  output  for  a  166 
Hz  tone  projected  from  the  source  described  above,  received 
on  the  single-line  HLA  array.  Beamformer  output  is  shown 
at  the  source  bearing  (endfire)  and  depth  (50  m),  and  is  nor¬ 
malized  by  a  noise  floor  estimated  from  the  25th  percentile 
of  all  beams.  With  no  peak-picking,  significant  range  am¬ 
biguities  are  seen  parallel  to  the  source  track.  These  ambi¬ 
guities  are  reduced  using  cluster-based  peak-picking.  If  a 
small  cluster  radius  is  used,  as  in  Figure  9,  the  rejection  of 
ambiguities  is  excellent.  However,  SNR  and  environmental 
mismatch  may  require  that  a  larger  cluster  size  be  chosen. 
Simulation  examples  have  shown  that  noticeable  cleanup  is 
also  obtained  for  a  cluster  radius  of  1.5  dB. 


Figure  9:  Simulated  results  for  cluster  peak-picking  with 
0.3  dB  radius,  single-line  HLA,  166  Hz  tone. 
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7.  SUMMARY 

Applying  matched  field  processing  concepts  to  sonar  arrays 
with  limited  aperture  often  results  in  a  set  of  highly  over¬ 
lapped  replica  vectors.  Clustering  methods  can  be  used  to 
identify  and  map  out  the  resulting  ambiguities  in  the  array 
manifold.  Ways  of  exploiting  the  cluster  map  have  been 
described  above.  As  demonstrated  using  SWellEx-96  HLA 
data,  bearings-only  clustering  can  show  mismatch  reduction 
benefits  in  comparison  to  direct-path  beamforming. 
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Matched  Field  Processing 
(MFP)  seeks  to  exploit 
shallow  water  multipath  for: 

-  3-D  source  localization 

-  Possible  detection  gains 

Leads  to  a  sparse  array 
problem  with  highly 
redundant  replicas 


Degrees  of  Freedom 
for  MFP  Replicas 
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